Friction in nanoelectromechanical systems: Clamping loss in the GHz regime 
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The performance of a wide variety of ultra-sensitive devices employing nanoelectromechanical 
resonators is determined by their mechanical quality factor, yet energy dissipation in these systems 
remains poorly understood. Here we develop a comprehensive theory of friction in high frequency 
resonators caused by the radiation of elastic energy into the support substrate, referred to as clamp- 
ing loss. The elastic radiation rate is found to be a strong increasing function of resonator frequency, 
and we argue that this mechanism will play an important role in future microwave-frequency devices. 

PACS numbers: 85.85.+j, 63.22.+m 



Nanoelectromechanical systems (NEMS) have enor- 
mous potential for both fundamental research and novel 
device applications 0, E H [1] • High-frequency NEMS 
resonators are the leading candidates in a current effort 
to observe quantum behavior in a macroscopic mechani- 
cal system Achieving the quantum limit requires 
resonator frequencies to that arc significantly higher than 
k^T/h, as well as high mechanical quality factors 



(a) 



Q = <-UT, 



(1) 



where r the energy relaxation time. In this limit the res- 
onator's phonons will be analogous to photons in an elec- 
tromagnetic cavity, and a variety of mechanical quantum 
"optics" effects have been discussed in the literature , in - 
cluding lasing j9(, squeezing 

Gnum, cooimg mni, 

and quantum nondemolition measurement [l5j. Macro- 
scopic quantum tunneling of a low-temperature beam 
has been predicted 0, 0], and methods to entangle 
NEMS with other quantum systems have been proposed 
13 El l^olj w hich may lead to their use in quantum 
computing IH S . NEMS operating in the 
classical regime have already been used successfully in 
several novel devices, particularly for sensing and actua- 
tion. They have been used to study the Casimir force at 
short distance s EilEp , and can be used as ultra-sensitive 
mass sensors |28l l29l|. recently detecting mass changes 
with zeptogram sensitivity [3(3 ] . They have demonstrated 
enormous potential in scanning probe microscopy [HI , in- 
cluding the ability to detect the spin of a single electron 
|.32j . and may also find future applications in classical 
signal processing Hi IM |M [H3| . 

The performance of a mechanical resonator improves 
with increasing Q. However, in a series of experi- 
ments on small mechanical resonators, performed over 
many years by several different groups, the quality fac- 
tor for motion in the fundamental flexural and other 
low-lying vibrational modes has been observed to de- 
crease approxima tely linearly with decreasing resonator 

size |aM EM EH, E3, [43, [43, E^, [43, E^. 

This dependence is usually interpreted to be a signa- 
ture of surface contamination effects, such as dissipation 





(c) 




FIG. 1: (Color online) Simulated resonators, (a) 2.4 GHz 
perpendicular beam resonator, (b) 3.2 GHz cantilever, (c) 
16 GHz doubly clamped beam. The material is Si and I = 
50 nm in the examples shown. The surfaces are colored ac- 
cording to the z component of their traction, with yellow rep- 
resenting zero. The lower surfaces are attached to a semi- 
infinite elastic substrate. 



caused by adsorbates or surface reconstructions. How- 
ever, measurements of the temperature dependence of Q 
also su gge st that glass-like structural two-level systems 
[finisHlfi^ lfiil or other internal "bulk" lattice [13 or elec- 
tronic |4fJ, |48| defects may also play an important role. 
Many energy dissipation mechanisms have been proposed 
to explain these experiments. However, most explana- 
tions are qualitative in nature, and a detailed quantita- 
tive analysis of the competing mechanisms has not yet 
been carried out. 

An obvious mechanism for energy dissipation in a me- 
chanical resonator is elastic radiation loss through the 
resonator supports, where the resonator is attached to a 
bulk substrate Any mechanical interaction 

of the resonator with its surroundings will allow vibra- 
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tional energy to escape to the substrate. But little is 
understood quantitatively about this so-called clamping 
loss mechanism, and some existing estimates are in con- 
tradiction with experiment. For example, the estimate 

M 



Qcl « 0.3 



(2) 



for an in-plane flexural mode applied to the 1 GHz SiC 
beam of Huang et al. results in a Q factor of 

about 230. smaller than observed. In the experiment 
of Rcf. [Sjj, the beam length is 1 = 1.1 /im and width is 
w=120nm 

In this paper we develop a general method for the accu- 
rate calculation of clamping loss rates for a wide variety 
of NEMS geometries. We assume at the outset that the 
clamping loss is small (its contribution to the total Q _1 
much less than unity), allowing a pcrturbative treatment. 
If the total dissipation rate Q^ 1 is also <C 1, as in the 
case of experimental interest, it will be given by a sum 
of contributions from each dissipation channel present. 
Our results suggest that clamping loss will play an im- 
portant role in future GHz-frequency devices, especially 
above 10 GHz. 

We begin our analysis with a discussion of the funda- 
mental mode |5flj of the perpendicular beam resonator 
illustrated in Fig. with frequency ui. The lower 
surface is to be attached to a semi-infinite elastic sub- 
strate, and all exposed surfaces are stress free. u> will 
always be degenerate with vibrational modes of the sub- 
strate, so how can clamping loss ever be small? The 
answer is that an elastic wave present in the beam will 
be strongly reflected at the substrate when its emitted 
wavelength lu/v, which for the fundamental mode is of 
the order of I, is larger than the transverse dimensions of 
the beam. Here v is a characteristic sound velocity of the 
substrate. The frequency-dependent transmission proba- 
bility between a semi-infinite wire (not finite beam) and 
a three-dimensional bulk solid similarly vanishes in the 
low frequency limit however, that transmission 

probability cannot be directly used to estimate clamping 
loss because it would only be valid for beam lengths large 
compared to the wavelength (to be able to neglect the ef- 
fect of cutting the wire to a finite length) , and hence only 
to the higher frequency non- fundamental modes. We cal- 
culate the clamping Q by regarding the resonator's vibra- 
tional eigenfunction stress distribution (calculated with 
a fixed lower boundary) to be a source of elastic radia- 
tion into the substrate. The radiation problem is treated 



exactly within continuum elasticity theory, and the res- 
onator vibrational modes arc calculated numerically us- 
ing finite-clement methods. Our method applies to a 
wide variety of (but not necessarily all) NEMS resonator 
geometries, and can be used at high fundamental-mode 
frequencies where the subdominant terms in a multipolc 
expansion would become important. 

The substrate is modeled as an isotropic elastic con- 
tinuum, which is sufficient for the frequencies of interest. 
The substrate occupies z > 0, and at the surface z = 
the traction is t(cc,y), which vanishes outside of the res- 
onator contact region. The elastic displacement field u is 
decomposed into longitudinal parts and transverse parts 
ui and u t , where 



ui(r) = 



d 2 K 
(2tt)2 



Ui(K) e iKr e iqiz 



and 



Here 



91 = \ / ~2 ~~ A ' 2 aIld <h = \ ~2~ K ' 



(3) 



(4) 



(5) 



K 

and Ut satisfy (K+q\e z ) xUi = and (K+g t e z )-U t = 0, 
and therefore can be written as Ui = (K+q\e z ) </>(K) and 
U t = 3>(K) — q^K. • 3>(K)e z . Here is a scalar field, 
and $ = & x e x + <S> y e y is a vector field. 

By construction, ui and u t satisfy the vector Helmholtz 
equations, conditions of longitudinality and transversal- 
ity, and have only outwardly radiating components (for 
real q). The elastic "potentials" cf> and $ are now fully 
determined by the three components of the surface trac- 
tion 



U = A(V • u)S iz + 2/iu lz , 



(6) 



with Uij the strain tensor and A, /i the Lame coefficients. 
In terms of the Fourier transform t(K) in the xy plane 
we obtain 



(7) 




where 
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S = 



[Kl{q% - Kl) - (q 2 - K*f - 4K 2 q m ]/q t K x K y (K 2 - q 2 + 4q m )/q t 2K M ' 

K x K y (K 2 - q 2 + Aq m )/q t [K*(q* - K 2 X ) - (q 2 - K 2 ) 2 - 4K 2 qi q t }/q t 2K y q iqt 

-2K x q t -2K y q t K 2 - q 2 ti 



(8) 



Hero x = n\{K" — It) + ^qiqtK ]. The clastic energy current, averaged over time, is 



j z = — Im (A u* z V • u + /i u* • Vu z + u* • d z u), 



(9) 



and the radiated power is 



P 



d rj z 



to 



pLJ' 



d 2 K 



{27t)2 |0| ffl + pt; t 



d 2 K 



t (2n) 2 



St * 



K ■ * 



St 



(10) 



TABLE I: Resonator frequency ui/2n and Q factor as a func- 
tion of size I. Beams have a 5nm thickness and 10 nm width. 





20 nm 


50 nm 


100 nm 


perpendicular 


14.7 GHz 


2.42 GHz 


607 MHz 


beam 


Q = 4.1xl0 3 


Q = 2.5xl0 5 


Q = 7.1xl0 6 


cantilever 


30.1GHz 


3.20 GHz 


695 MHz 


beam 


Q = 1.0xl0 3 


Q = 4.1xl0 5 


Q = 1.4xl0 7 


doubly clamped 


97.1GHz 


16.0 GHz 


3.85 GHz 


beam 


Q = 21 


Q = 2.1xl0 3 


Q = 7.9xl0 4 



The integration domains T\ and I\ in (|l()|l have |K| less 
than uj/vi and tu/vt, respectively. Given the surface trac- 
tion t(x,y), P is now readily evaluated. The damping 
time r is E/P, with E the initial energy of the resonator. 
The required traction information can be accurately de- 
termined by numerical finite-clement methods, without 
the need to include a large "environment" in the simu- 
lation. Our numerical results converge to the precision 
reported here after a few tens of thousands of volume 
elements. 

Our results are summarized in Table [I] The Q factor 
for a perpendicular beam resonator shows that clamp- 
ing loss in that case is probably negligible below 1 GHz, 
but that it becomes quite significant above 10 GHz. An 
approximate analytic expression for Q in this geometry 
can be obtained by constructing the fundamental-mode 
cigenfunction with thin-beam theory, and keeping only 
the leading (monopole) transverse stress moment in l)l()|l. 
This leads to 

where (3 is a material-dependent parameter, equal to 
0.194 for Si, and where 

Q = izv\A~ 1/2 , with A = tw. (12) 



TABLE II: Asymptotic Q factor 1111 for the perpendicular 
beam resonator. Dimensions are the same as in Table [I] The 
frequency is calculated analytically using thin-beam theory, 
and n/27r is 599 GHz. Q PJ is taken from Ref. |H|]. 



perpendicular beam 


20 nm 


50 nm 


100 nm 


frequency 


17.1 GHz 


2.73 GHz 


683 MHz 


asymptotic Q 


2.5 xlO 3 


2.2 xlO 5 


7.1 xlO 6 


Qpj 


1.7xl0 3 


1.6 xlO 5 


5.2 xlO 6 



fl is a characteristic frequency above which the resonator 
beam — with cross-sectional area A — can support nonuni- 
form transverse modes, changing the dynamics from one 
dimensional to three dimensional. The validity of the re- 
sult (jl 1|) relies on the applicability of thin-beam theory; 
it is therefore asymptotically exact in the u> Q limit. 
The beam-theory frequency and asymptotic Q factor for 
the perpendicular resonator is given in Table HTl showing 
that the expression l|ll(l becomes unreliable (overestimat- 
ing dissipation) above about 10 GHz, which is the regime 
where clamping loss just begins to be important. Also 
given in Table ^ is a simpler (but not asymptotically 
exact) estimate by Photiadis and Judge [5rj| . 



which predicts clamping loss rates somewhat higher than 
that given by (fTTf l . 

The Q factor for the cantilever and doubly clamped 
beam resonators behave quite similarly to the perpendic- 
ular resonator, when viewed as a function of frequency 
instead of size. In all cases the Q factors at a few GHz arc 
about 10 5 , but fall to around 10 3 in the few tens of GHz 
regime. Loss at 100 GHz is severe. It is reasonable to 
assume that this trend is quite general, and will apply to 
other resonator geometries as well. The development of 
ultrahigh-Q NEMS resonators above about 10 GHz will 
require phononic bandgap engineering |6ll| or other meth- 
ods of acoustic isolation. 
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Finally, we briefly comment on the high frequency be- 
havior of the asymptotic expression (|llf> . which we have 
argued to be exact (for the perpendicular beam res- 
onator) in the low-frequency limit. It follows from our 
earlier discussion that Q would be expected to be of order 
unity at the frequency fl Thus, it is tempting to conclude 
that the formula can be used for all frequencies in 
the range < u) < i7, which would be incorrect, despite 
the fact that it is accurate at the two endpoints. 

This work was supported by the NSF under grants 
DMR-0093217 and CMS-040403. 
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